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The decay rate of metastable states is determined at high temperatures by thermal activation, 
j^jfj' whereas at temperatures close to zero quantum tunneling is relevant. At some temperature T c 

r-{ i the transition from classical to quantum-dominated decay occurs. The transition can be first-order 

like, with a discontinuous first derivative of the Euclidean action, or smooth with only a second 
derivative developing a jump. In the former case the crossover temperature T c cannot be calculated 
perturbatively and must be found as the intersection point of the Euclidean actions calculated at 
low and high temperatures. In this paper we present a sufficient criterion for a first-order transition 
in tunneling problems and apply it to the problem of the tunneling of strings. It is shown that the 
. problem of the depinning of a massive string from a linear defect in the presence of an arbitrarily 

D ' strong dissipation exhibits a first-order transition. 
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a Investigations of the decay rate of metastable states have a long history, going back to KramersS who calculated 

the lifetime of a classical trapped particle separated from the true equilibrium state by a high potential barrier n Since 
then different decay phenomena have been investigated. The motion of dislocations across the Peierls barrier J3tl the 
decay of the current in a Josephson loop,cl and the creep of vortices in superconductorscl are typical and well-known 
examples of metastability. At sufficiently high temperatures a metastable system decays due to thermal activation and 
the decay rate T obeys the Arrhcnius law T ~ exp (-U/T), whereas at temperatures close to absolute zero quantum 
tunneling is relevant and T ~ exp (-S/K), with 5* the action at zero temperature. Below we shall consider systems 
for which the semiclassical description is applicable. In this case U/T and S/Ti are large; otherwise, the system would 
not be truly metastable. 

Let u(t, r) denote the coordinates of the system under consideration depending in general on imaginary time t 
, and spatial variables r. Within the semiclassical approximation the decay rate at a temperature T is determined by 
■ the contribution of the trajectories close to the one extremizing the Euclidean action and satisfying the periodicity 
t^- I condition u(0, r) = u(h/T, r). In the high-temperature regime the function extremizing the action is time independent 
ON ■ and, consequently, the activation barrier in the Arrhenius law U does not depend on the dynamic properties of the 
system, meanwhile the bounce solution at low temperatures is time-dependent and the dynamics enters. At some 
temperature T c the transition from the time-independent to the time-dependent solution occurs. In some cases the 
bounce solution just below the crossover point can be written in the form 

I 

u(i,r) = u t/l (r) + <Su(r)cos^^^J, (1) 



I. INTRODUCTION 



where the function <5u(r) is considered to be small. Such an assumption leads to an action Se uc \(T) with a continuous 
first derivative in the crossover point (see Fig. 1). 

The second derivative of S is discontinuous at the point T c . The second derivative has a jump at T = T c . Following 
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Ref.EI we shall call this situation a "second-order transition in the crossover point." In this case the bounce solution 
with the minimal Euclidean action at low temperatures can be obtained by a continuous deformation of the thermal 
solution. However, in general such a deformation is not possible: The trajectory corresponding to the minimal action 
may jump at a certain temperature. In this case we deal with a first-order transition: it can be shown that the first 
derivative of S-e uc \{T) has a jump at the transition point (see Fig. 2) and the expansion (|l|) is not valid. 

Strictly speaking, even if we deal with a first-order transition in a tunneling problem, there is a narrow crossover 
region from one solution to another because contributions originate from several saddle points and we need to take 
into account all of them. However, the better the semiclassical approximation is applicable, the narrower this region 
becomes. 
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The problem considered here is related to the mean-field theory of phase transitions. The Euclidean action can be 
identified with the free energy in the Landau theory with an order parameter 5(r, T) defined as 5{r, T) — max t |u(t, r) — 
Uth(r)|, for example. At high temperatures u(i, r) = u t / l (r) and S(r,T) = 0, whereas at T < T c , 5(r,T) ^ 0. The 
order parameter changes continuously at the point T = T c in the case of a second-order transition and discontinuously 
if a first-order transition takes place. 
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Fig.l Euclidean action S as a function of temperature T for the case of a second-order transition. 
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Fig. 2 Euclidean action S as a function of temperature T for the case of a first-order transition. The first derivative of S is discontinuous 
at the point T c . 



In this paper we present a simple criterion for the appearance of a first-order transition from quantum to classical 
behavior in tunneling problems and apply it to the problem of the tunneling of strings. 



II. OSCILLATIONS CLOSE TO THE THERMAL SADDLE-POINT SOLUTION 



In this section we study the behavior of the imaginary time oscillation period of the solution of the equation of motion 



close to the thermal saddle-point solution. In section II A we describe the 



sneral theory for ape-dimensional (ID) 
] and Chudnovskytl for ID Hamiltonian 



Hamiltonian systems. We briefly summarize the results of Lifshitz and Kagar 
systems and provide an expression fo r the jump of the derivative of the Euclidean action (see Eq. (|13|)) for the case 
of a first-order transition. In section II B we use a perturbative approach for the calculation of the imaginary time 
oscillation period in the vicinity of the thermal saddle-point solution, and define a criterion for a first-order transition. 



A. General theory for ID Hamiltonian systems 

As shown by Chudnovskyi, if the imaginary time oscillation period of a massive particle is not a monotonous 
function of energy, a first-order transition from quantum to classical behavior takes place. This statement remains 
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true for any one-dimensional Hamiltonian system: Consider a ID system subject to potential which allows tunneling 
in one direction (Fig. 3). 




u u 

Fig. 3 ID potential well with a metastable state. A particle trapped at u = can tunnel out to the right in order to reach the true 
ground state. 

We choose the metastable minimum of the potential as our zero energy. Let us define the height of the potential barrier 
as Uq. Using the semiclassical approximation we can write the imaginary parts hT n of the energy E n of the eigenstates 
in the forma 



hT r . 



fioj{E n ) 



4tt 



■ exp 



(2) 



where p is the generalized momentum, E n < Uq are the energy levels in the equivalent potential without tunneling, 
and Lo(E n ) are the corresponding classical oscillation frequencies. At finite temperature T the decay rate can be easily 
found by averaging over the Boltzmann distribution 



2 x - uj(E„) ( E n S 1 ., 



4tt 



T 



(3) 



where Z = ex P (~ ^f-) is the partition function of the particle in the well. With Eq. (||) we reproduce the 
well-known expression Y = (2/h)lmF for the decay rate at low T (see RefJlj). 

Equation (||) is applicable only at sufficiently low temperatures. At high temperatures the excitations with energies 
larger than the barrier .height are relevant in the determination of the preexponential factor. The decay rate then is 
given by the equationO'EJO 



r 



uj(Eq) 
2tt 



(4) 



where uj(Eq) is the oscillation frequency of the system near the metastable minimum, ft is necessary to emphasize 
that the above expressions for the decay rate are applicable in the case of intermediate friction r\. E.g., for a particle 
of mass m the criterion ujT/lL.. <ti w /m <C u>, where u> 2 = —U"(uo)/m with uq the position of the maximum of the 
metastable potential (see Refs.ElEilii3), guarantees, that on the one hand the system in the well is properly equilibrated, 
while on the other hand the dynamics is not affected by the damping. 

The sum in Eq. (||) can be well approximated by the term where the function fr(E) = —E/T — S(E)/h takes its 
maximal value and the problem reduces to the calculation of the maximum of Jt{E) within the interval [0, Uq]. The 
extremal condition for fr(E) takes the form 



dE 



-1-1^ = 0. 
T hdE 



(5) 



It is well known from standard courses of classical mechanicsEl that for a periodic problem the derivative of the 
action with respect to the energy is equal to the oscillation time r at this energy. In our case the energy axis has the 
opposite direction (the action S(E) corresponds to the motion in the inverse potential with energy E). Hence, 



^- = -l = -r(E), 
dE T V h 



(6) 
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i.e., the decay is dominated by the Euclidean action with the trajectory periodic in the time interval h/T. 

If the function t(E) is a monotonously decreasing function of energy, Eq. (^|) has one solution for T < Ti/tq 
(tq = t(Uo)) and no solutions for T > Ti/tq. At temperature T c = Ti/tq a second-order transition takes place. 

Next, let us suppose that t(E) goes through a minimum in the interval (0, Uq) (see Fig. 4). In this case it can be 
showntl that there is an energy E\ € (0, E m [ n ) : and associated with it a temperature T c = Ti/t(Ei), where a first-order 
transition takes place. 

Let us show that if E\ ^ Uo, the Euclidean action Seuci(T) determining the decay rate T indeed has a discontinuous 
derivative at the point T — T c (see Eq. ([§). At T > T c 

^thermal 

( T ) _ U o 

n n t ■ u 




Fig. 4 Non-monotonous dependence of the oscillation time r in the inverse potential on the energy E. The arrow indicates the jump of 
the semiclassical trajectory from the thermal assisted quantum solution to the classical one. 



(8) 



The action just below T c can be written as (see Eq. (|3|) 

SEucl(T) ^quantum (T) _ E X {T) S(E{) 

h h T n 

(note that E\ is a function of temperature). 

The variations of ^thermal and ^quantum with respect to the T near the point T c are 

crmal 

Z = , 

h T r 2 



(9) 
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Ei(T c ) fj_,ldS_ 
T r 2 \T C + HdE 



E=E 1 (T C ) 



V-^ST. 
J dT 



Taking into account that at the point E\ 

dfr 
dE 



E=E X 



1 



1 dS 
hdE 



E=E X 



we obtain 
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quantum 



T r 2 



ST. 



Thus, the jump of the derivative at the point T c is 

dS-Eud h (Ei - U ) 



A- 



dT 



T c 2 



<0, 



(10) 



(11) 



(12) 



(13) 



and the thermal action always decays more rapidly than the quantum one. We thus have shown that if the imaginary 
time oscillation time r is not a monotonous function of energy E, a first-order transition takes place. Equation ( p^ ) 
quantifies the strength of the first-order transition: if E\(T C ) ~ Uq, the difference in the bounce solutions is small and 
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the first-order transition is weak. In the opposite with case Ei(T c ) ~ the saddle-point solution deforms significantly 
and we have a strong first-order transition. 

For a ID massive particle the above physics can be easily realized. Consider a given function t(E) for a trapped par- 
ticle. Then there exist infinitely many potentials reproducing t(E) as their oscillation timesEj. A simple way to obtain 
a non-monotonous behavior is to choose a function t(E) with one minimum. For more complicated nonmonotonous 
dependences the system can exhibit several transitions. In general a system may exhibit several first-order transitions 
but not more than one of second-order. 

All the results obtained in this section remain true for an arbitrary ID metastable Hamiltonian system as Eq. @ 
is applicable in this case. 



B. Nonlinear perturbation near the thermal saddle-point solution 

Consider a metastable system whose oscillation period in the imaginary time r is a function of a scalar parameter a, 
t = r(a). For Hamiltonian systems it is convenient to choose this parameter to be equal to the energy, i.e., a = E, see 



section [I A. However, for dissipative systems one cannot use the energy as it is no more a conserved quantity. In this 
case we can parametrize the periodic imaginary time solutions of the equation of motion by the amplitude a of the 
oscillations: The amplitude a quantifies the difference between the thermal and the weakly time-dependent solutions 
in the vicinity of the thermal saddle-point solution (later we will show how to define such a notion as the "amplitude" 
for any metastable system). It can be easily understood that if r(o) is not a monotonous function of the amplitude, 
the system exhibits a first-order transition: Starting from the zero temperature bounce solution the period r = h/T 
decreases with increasing temperature. However we cannot carry this solution beyond the temperature h/T m i n and 
thus will encounter a first-order jump to the thermal solution at some temperature T c < h/T min (in the simplest case 
of a r-dependence with one minimum). 

The basic idea leading to the criterion for a sharp transition from quantum to classical behavior is the following: 
Let us investigate the imaginary time oscillation period in the vicinity of the thermal saddle-point, where a = 0. From 
the above discussion follows that if r(a) — r(0) < 0, the system will exhibit a first-order transition. It turns out that 
a general expression for r(o) can be obtained for a wide class of systems. 

At high temperatures the saddle-point solution is time-independent. Let us slightly perturb this solution and 
calculate its oscillation time in the vicinity of the point a = 0. We suppose that the imaginary time Lagrangian L of 
the system under consideration can be written in the form 

L = T(u,u,r) + V(u,r), (14) 

where the term T(u, u, r) is responsible for the dynamical properties of the system and V(u, r) is the potential energy. 
The equation of motion corresponding to the Lagrangian ( pL4|) takes the form 

fu=£. (15) 

Later we shall suppose that the operator I is linear and that Sl/6r — 0. This is not a strong restriction as the dynamical 
terms traditionally considered in the Lagrangian (fhuLconsist of dissipative, massive, and Hall terms. Dissipation is 
usually described by the Caldeira-Leggett formalisms, leading to a linear term in the equation of motion. Massive 
and Hall terms also lead to linear equations, i.e., the operator I satisfies the conditions discussed above. For a massive 
particle the operator I takes the form I = md 2 /dt 2 . 



In the high-temperature regime the solution of Eq. ( TL5h u t ^(r) is time-independent. Let us expand Eq.(15) into a 
series around this solution. Substituting u(r, t) — Uth(rJ + Su(r, t) into Eq. ( J15| ) and expanding in <5u we obtain 

ISu = hSu + G 2 [Su] + G 3 [Su], (16) 

where I = S 2 V/6u 2 is a linear operator and G 2 and G3 arc nonlinear operators satisfying the conditions 

G 2 [Ay] = A 2 G 2 [y] and G 3 [Ay] = A 3 G 3 [y], (17) 

where y is an arbitrary vector and A is a constant. Our goal is to solve Eq. ( [l6| ) for Su(r, t) and find the correction 
to the oscillation period away from the thermal saddle point. To lowest order in perturbation theory we separate 
variables with the Ansatz 

Su = au (r) cos uj^t (18) 
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and substituting (|18[) into (|16|) while neglecting terms of order higher than a and higher we obtain 



u o(r) 



I COSLUot 



hu (r) 



COS U>(jt. 



(19) 



The function coswo^ is an eigenfunction of the operator I, Icostoot = I (loq) cos Lo^t, and we obtain the eigenvalue 
equation for u (r) 



l(u )u (r) = hu (r). 



(20) 



The operator h has one negative eigenvalue due to the unstable direction in phase space. On the other hand, 
1(loq) < 0, i.e., Eq. ( pp| ) has only one solution. The function uo(r) then is the eigenfunction of the operator h 
corresponding to its lowest eigenvalue ho and the frequency loq is determined by the condition l(u>o) — ho. Note that 
T c = hLOo/2ir is the transition temperature in case of a second-order transition. At this stage it becomes clear how 
to define such a notion as an "amplitude" for any tunneling problem: it is the expansion coefficient of the unstable 
mode. 

Next, let us write 



5u = au (r) cosuit + Su 2 (r, t) , to = io + ^2, 



(21) 



where w 2 is the correction to the frequency loq and 5u 2 ~ a 2 . Substituting (ElJ) into (16), neglecting terms of order 



and higher, and using ( |l7|) we obtain the equation for <5u 2 

Z(w)au (r) cos tot + I6u 2 = /(wo) au o( r ) coscut + hSu 2 + a 2 G2 [u (r)]cos 2 cut. 
Rearranging terms, we arrive at 



5u 2 



(Z(wo) - auo(r) cos (u>t) + — G 2 [u (r)] (1 + cos (2wt)) 



(22) 



(23) 



Since (I — /i)uo(r) cos Lot = 0, the first term has to vanish and we obtain no shift in the frequency, l(u>) = 1(loo), i.e., 
u> = ujq, L02 = 0. The solution of Eq.(E2n then reads 



where 



^u 2 = gi (r) + g 2 (r) cos 2uj t , 



gl (r) = --a 2 h- 1 G 2 [u a (r)}, 



(24) 



(25) 



g 2 (r) = - l -a 2 {h-l{2u Jo )) 1 G 2 [u (r)}. 

To third order in perturbation theory we make the Ansatz 

Su = auo(r) cosut + <Su 2 (r, t) + #113 (r, t), to = u>o + 0)3. 

Substituting (^7| ) into ( |l6| ) and neglecting terms of order a 4 and higher we obtain the equation for 6113, 

/(w)auo(r) cosiot + Z(2cj)g 2 (r) cos 2iot + ISu^ 
= Z(wo)auo(r) cos tot + hgi + hg 2 cos2iot + h5u^ 
+ G 2 [auo(r) cos Lot + gi + g 2 cos2tot] + G3 [auo(r) cos tot] . 



(26) 



(27) 



(28) 



As gi (r) , g 2 (r) ~ a 2 , we can expand the term G 2 [au (r) cos tot + gi + g 2 cos 2iot] around the function au (r) cos tot : 



G 2 [au (r ) cos tot + gi (r) + g 2 (r) cos 2Lot] 



SG 

G2 [auo (r ) cos Lot] + a cos tot— — 

du 



(gi(r) +g 2 (r) cos2wt) 



(29) 



u=u (r) 
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Note that the opeator <5G 2 /<5u| UQ is a linear operator, see Eq. (|i~7j). 

As before we sum up all the resonant terms ~ coswiu (r) and obtain the following equation for the shift in the 
oscillation frequency 

1(cj) - l(u ) = (u , f Mr)]) , (30) 

(,u , u J 

where 

G 2 [u (r)] + ^G 3 [u (r)] (31) 

and (x, y) denotes a scalar product (see Appendix for a general discussion) . The criterion for a first-order transition 
reads r(a) — r(0) < and using dl/dui < we find 

(uo(r),f[uo(r)])<0. (32) 

In the next section we apply the criterion (|32] ) to the tunneling of strings. 



1 <5G 2 

f[uo(r)] = -2^7 



In (r~) 



III. TUNNELING OF STRINGS 



A. General theory 



In this section we apply Eqs. (|3l| ) and (|32] ) to the problem of a driven string i) tunneling between two potential 
wells and ii) depinning from a columnar defect. Within the linear elasticity theory the imaginary time Lagrangian of 
the string takes the form 



L/2 



c 



dz lL D (u,d t u) + U^ ) +T(V)- Fu 



(33) 



-L/2 



Below we consider the dynamical term L£>(u,d t u) to be the sum of massive and dissipative terms (we ignore the 
transverse contribution from the Hall term) 



+H/2T 



L D (u, d t u) = ^{d t uf - ^-d t u 



dt\ In 



-h/2T 



T 



dt,u. 



(34) 



L and e are, respectively, the length and the elasticity of the string, V(u) is the potential, the external force is assumed 
to be small: F <C F c with F c the depinning force. The thermal saddle-point solution Uth{^) corresponding to the 
Lagrangian satisfies the equation 



d 2 u _dV 
dz 2 du 



(35) 



The operators h, G2, G3, and 6G2/8u\ Uo ( z ) take the following form, 

f d 2 d 2 v 



dz 2 du 2 



(36) 



Uth(z) 



1 d 3 V 



2 du 3 



v\ G 3 [y] = - 



ld 4 V 



6 du 4 



SG 2 
Su 



(y) 



d 3 V 



du 3 



uo(z)y, 



(37) 



where y is an arbitrary function of z. The scalar product (2/1,2/2) of two arbitrary functions yi(z) and 2/2(2) is defined 

aS I-L/2 2/l(^)3/2 {z)dz. 



L/2 

Eq. (|32|) takes the form 
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L/2 

1 f d 3 V 
(u ,/N) = -| / dz-^ 



-L/2 



L/2 

1 f 1 d 4 V 

dz 



du 4 



ul < 0, 



-L/2 



-i\ d 3 V 



du 3 



(38) 



where l{w) = —puj 2 — ri\u!\. It is appropriate to mention that the operator h has one zero eigenvalue in the limit 
L — > oo, (sec III B| ) such that taking its inverse hr 1 needs some care: The eigenfunction corresponding to the zero 
eigenvalue is an odd function of z, whereas u 2 ] d 3 V / du 3 \ u is an even function and there is no contribution to the 
integral in Eq. ( |38| ) . Consequently, the contribution of the eigenfunction corresponding to the zero eigenvalue is equal 
to zero. 

Let us investigate the sign of Eq. (138) in the two cases mentioned above. 



B. Motion of a string across a slightly tilted periodic potential 

The tunneling of a string across a titled periodic potential is a convenient model for the description of the motion 
of dislocations across the Peierls barrier. In the case being studied V(u) is the periodic part of the potential (see 
Eq. (H^)). Below we denote time by t and the oscillation period by r. This problem shows a smooth crossover from 
a purely thermal to a thermal-assisted quantum activation at a temperature T c = T of order h(F/ pd) 1 ^ 2 (see Ref.i) 
with d the period of the potential V(u). The schematic evolution of the nucleus is shown in Figs. 5(a)-5(c). At 

1 /2 

temperatures smaller than T\ ~ UF/ (pV ) , with V the amplitude of the potenrial, the nucleus is of circular shape 
(see Fig. 5(a)) 






(a) 



(b) 



(c) 



Fig. 5 Evolution of the nucleus (in coordinates z = z/y/e and t = t/yTp). (a) Nucleus at low temperatures, (b) Nucleus at intermediate 
temperatures where the interaction between the nucleus' walls is relevant, (c) Nucleus at high temperatures. The string overcomes the 
barrier due to thermal activation. 



with the bounce solution u(z, t) = aO(R — r), R ~ y/Va/F is the radius of the nucleus (r 2 = z 2 /e + t 2 / p = S 2 + t 2 ) and 
the Euclidean action S ~ ^/pedVo/F. At T > T\ we have to take into account the interaction of the nucleus' walls 
(Fig. 5(b)). At T > T c — T the string overcomes the barrier due to thermal activation with U ~ d^/eVo (Fig. 5(c)). 
The thermal solution then can be obtained by a continuous deformation of the zero temperature bounce. As the 
problem exhibits a second-order transition from quantum to classical behavior, we should expect 1(uj) — l(u>o) > 0. 
Let us prove this inequality. The "potential" in the one-dimensional Schrodinger operator h (see (|36|)) is shown in 
Fig. 6. The spectrum of h consists of one negative, one zero (in the limit L — > oo), and positive eigenvalues (one can 
easily verify that the function duth/dz is an eigenfunction of the operator h with zero eigenvalue; on the other hand, 
duth/dz has one node, i.e., there is one eigenfunction with a negative eigenvalue). In the limit F —> the negative 
eigenvalue ho tends to zero as —F, see Ref.cl. The positive eigenvalues of the discrete spectrum of the operator h are 

of order Vo/d 2 . Analysing the terms in Eq. (|38|) we note that the operator 



1 



/2(h-l(2w ) 



has a large 
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negative eigenvalue ~ — 1/F originating from the unstable direction in the phase space. The remaining spectrum is 
positive and does not depend on F, as is also the case for the other factors appearing in Eq. (|38|). 




Fig. 6 The potential in the ID Schrodinger operator for the linearized problem of a string tunneling of a string across a slightly titled 
periodic potential. Uu n (z) = d 2 U /dz 2 \ Uth ( z y The negative eigenvalue appears because of the interaction of the two wells. The distance 
between the negative and zero eigenvalues is much smaller than the distance between possible positive eigenvalues of the discrete spectrum 
and the zero-eigenvalue level. 

Thus we can approximate (itn, / [mo]) by the contribution from the unstable direction and arrive at the estimate 



{u ,f[u }) 



h 2{h -l(2w )) 



( L/2 



(39) 



\-L/2 



The second term in Eq. ( p8| ) does not involve the small parameter F and can always be neglected in the limit F — > 0. 
As 1(uj) = —puj 2 — rj\uj\, ho = 1(ujo) we obtain 



1 



h 2(/i -/(2w )) 



< 



(40) 



and, consequently, l(u>) ~ l(u>o) > 0. We have shown therefore that in the vicinity of the point a — 0, r(a) is an 
increasing function of a. This is in agreement with the fact that the problem of the tunneling of a massive string 
across a slightly titled periodic potential exhibits a second-order transitional. 



C. Depinning of a string from a linear defect 

The behavior described above should be contrasted with the sharp transition obtained recently for the tunneling of 
a massive string from a linear defectO. The imaginary time Lagrangian has the form ( |33| ) (j] = 0) with the potential 
U(u) describing a single potential well of depth Vo and radius d. At low temperatures the nucleus has a circular shape 
as well. The macroscopic part of the bounce solution can be written in the form 

u{Z)t ) = ~r 2 + ^, r 2 = - + - = z 2 + P. (41) 
4 t e p 



The substitution of Eq. (|4l|) into Eq. (p3|) gives the Euclidean action 



fL[u( z ,t)]dt=^(^) . (42) 



h / 1 r " J h \ f 



The thermal exponent is given by U/T = (4^2/3) ^/e\^(V / FT). S q /h and U/T become equal at T c = 

(V2/3Tr)hF / (pVo) 1 ^ 2 ■ On the other hand, the dynamical solution ( f4l| ) remains valid as long as the radius of the 
nucleus is smaller tha n the periodicity in t, i.e., for T < T\ = 3n/8T c . The formal application of the perturbative 



procedure (see section II B) gives the following linearized problem determining the crossover temperature T , 
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d 2 u d 2 V 



dz 2 du 2 



{2irT 

U0 = ^It 



u . 



(43) 



The ID Schrodinger operator in (E3) has_again one negative eigenvalue which does not depend on the details of the 
potential V(u) in the limit F -> 0(Ref.BEa). T then is givep-by the formula T = ( f i/2 3 / 2 n)hF j Jp% = 0.900T C 
with fj, = 1.199 the root of the equation /itanh/i = 1, see Ref£3. As Tq < T c < Ti, the string exhibits a first-order 
transition at the temperature T c . In this case the quantum solution (see Fig. 5(a)) jumps at T c to the thermal solution 
(Fig. 5(c)). The intermediate regime (Fig. 5(b)) is absent. 

This sharp transition is again in agreement with the above considerations: There exist solutions of the equation 
of motion in imaginary time with periods r g (h/Ti,oo) and energy close to zero (we consider a massive dynamics 
here). On the other hand, at the temperature Tq there appears a time-dependent solution of the equation of motion 
with a period h/T and an energy equal to the energy of the thermal solution U. As h/T\ < h/T , the dependence of 
the oscillation period on energy is not monotonous and a first-order transition takes place. In the present case the 
jump of the derivative [Ei(T c ) — U]/T 2 in the action is large: Ei(T c ) is close to zero, U and TjT 1 are proportional to 
1/F, hence A(<9Seuci/<9T) ~ — 1/F 3 and one can see that we deal with a strong first-order transition. 

Let us find the sign of l(u>) — l(u>o) in (30) for the case r\ ^ 0, p 7^ 0. The "potential" in the one-dimensional 
Schrodinger operator h is shown in Fig. 7. 



Vo 
d 2 



U.(z) 



~F 2 



~F 2 




D-1/F 



Fig. 7 The potential in the ID Schrodinger operator for the linearized problem of the dcpinning of a string from a linear defect. 
Uu n (z) = d 2 U /dz 2 \ Uth ( z y The modulus of the negative eigenvalue is of the same order as positive eigenvalues in the discrete spectrum. 



In the limit F — > the most significant contribution to the integrals in ( pq ) arises from the region \z\ > D = y/2eUo/ ' F. 
The spectrum of the operator h consists of one negative, one zero, and positive eigenvalues. In the limit F — > the 
negative eigenvalue and the lowest positive eigenvalues of the descrete spectrum are proportional to F 2 and we need 
to analyze the first integral in Eq. ( |38| ) more carefully than before. The second integral in Eq. (^8|) is again irrelevant. 
In the limit F — > 0, the eigenfunctions of the lowest levels take the same form in the region \z\ > D, i.e., we can write 



("o,/K]) ! 



L/2 



dz 



d 3 V 



-L/2 
L/2 



dz 



du 3 
d 3 V 



h - 1(2uj ) 



d 3 V 



du 3 



du 3 



-L/2 




^(h n -l{2u ))- 1 



\u n )(u r , 



d 3 V 



du 3 



dzu n 



d 3 V 



du 3 



(44) 



iz\>D 



Here, we have used that d 3 V/du 3 
normalizers. Consequently, 



for \z\ < 0. In the region \z\ > D, u n = ±C„Mo/Co with C n denoting the 
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(uqJ[uq\) = - 



d 3 V 



<9u 3 



\z\>D 



UqcIz 



(45) 



For the lowest eigenvalue we have uo(|z| < D) = Co cosh —^zj and the eigenfunction belonging to the zero 

eigenvalue takes the form Mi(|z| < D) = C\z. In the limit F — > we have u' n /u n \ z= o- = u' n /u n \ z =D+ = u 'a/ u o\z=D+ , 
i.e., u' n /u n \ z =D- = Uq/uo\ z= d- ( we use the condition of continuity of the logarithm derivative), i.e., for the lowest 
eigenvalue we obtain the condition 



-D tanh 




1. 



(46) 



Only even levels contribute to the integral in ([38|) and their eigenfunctions are given by C n cos ( y ) , hence we 
obtain for h„ 



D tan 




The most significant contribution to the integrals J_ 

1 



L/2 



L/2 



1. 



\dz arises from the interval \z\ < D, i.e. 



n 2 — 



D 

J cosh 2 

-D 



— z)dz 



D cosli 



and 



D 

J COS 2 

-D 



^z)dz -Dcos 2 



D 



2,4. 



(47) 



(48) 



(49) 



In Eqs. (|4J) and (|9|) we have used Eqs. (|46|) and (|7|). 

Let us suppose that r\ = 0. Later we will show that if the problem exhibits a first-order transition for r\ = 0, 

the kind of transition remains the same for any r\ ^ 0. Numerical solution of Eqs. ([l6|) and ( fi"7| ) show that ho = 

-1.439229e/L> 2 and h 2 = 7.830964e/L> 2 . For the case of a purely massive dynamics 1(2uj ) = AI{uj ) = Ah a . The 

-l 



(50) 



operator h + ( 1/2) f h — 1(2luq) ) has only one negative eigenvalue, i.e., if we show that the sum 



S 



1 



M 



1 



0,2.. 



2(fc n -i(2w )) < 

is positive for some finite M, a first-order transition takes place. With M = 2we obtain 



i 



i 



fto 2(/i -Z(2o;o)) 



1 



2(^-^(2^)) 



0.008796D/e > 0, 



i.e., the transition from quantum to classical behavior is of first-order. 

Finally, let us show that the kind of transition remains the same for any 77 ^ 0. We can write for the ratio 



1(2lu ) _ 4/Ocjo 2 + 2rjLU 
1(uj q ) 



puj 2 + rjojo 



= 4- 



2rjuj 



i.e., the relative contribution of the operator (1/2) (h — 1(2ujo 
of the operator h^ 1 . 



-1 



<4, 



(51) 



(52) 



becomes even larger compared to the contribution 



Comparing the problems considered in sections III B and [II C the main difference appears because of the different 
character of the spectrum of the linearized problem. For the case of a string tunneling across a slightly titled periodic 
potential the most significant contribution arises from the negative eigenvalue. The contributions of the positive 
eigenvalues are negligible in the limit F — > 0. On the other hand, for a string depinning from a linear defect the 
contributions of the positive eigenvalues are of the same order as that of the negative eigenvalue. 
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IV. CONCLUSION 



If a second order transition takes place, it is possible to use a perturbative approach for the calculation of the 
crossover temperature. In this case we have to substitute equation (|l|) into the equation of motion. Close to T c the 
function Su(r) is small and we obtain the linearized eigenvalue problem (see also Eq. (j2T 

£ 2 Seuci[u] 



5u 2 



8u = XSu, (53) 

u=uo(r) 



where the only negative eigenvalue of Eq. ( p3| ) (due to the single unstable direction in the phase space) determines 
the transition temperature T c . For the case of a first order transition at T c the formal application of the perturbative 
procedure again gives some "crossover" temperature Tq which comes to lie below the true transition at T c , however. 
For the dependence of the oscillation time on energy considered above (see Fig. 4) To — %/tq. In Fig. 8 we plot the 
actions for the various extremal solutions corresponding to the situation in Fig. 4. 




To T c Tmin T 

Fig. 8 Euclidean action S as a function of temperature T for a non-monotonous oscillation time t(E) as given in Fig. 4. 

The solid line AC corresponds to quantum tunneling. At the point C, a first-order transition takes place and the 
solid line CB shows the thermal action. In the point E the "false" crossover takes place. The line EF corresponds 
to the oscillating motion of the tunneling system with energy E £ (E m i n ,U~o) (see Fig. 4). The dashed line CF is 
obtained by the continuation of the quantum action into the region T > T c and corresponds to the motion with energy 
E E (Ei(T c ), E m i n ). In the point F this line and the "false" crossover action intersect_each other. The dashed line 
EC is a formal continuation of the thermal action into the region T <T C , see also RefJa. 

The above discussion shows that in the case of a first-order transition one cannot use the perturbative representation 
of the bounce solution (see equation (Q)). The transition takes place at a temperature higher than that given by the 
formal solution of the linearized eigenvalue problem. In this case one has to find the low temperature expression for 
the Euclidean action and determine the point where it is equal to the thermal action. 

There are a number of articles concerning the decay of metastable states, where a sharp transition from quantum, 
to classical behavior has been discovered. A first-order transition has been obtained by Lifshitz and Kagan in Ref.EJ 
concerning phase transitions at low temperatures, by Chudnovskytl for a massive particle in a ID potential for the 
case of a nonnmpotonous dependence of the imaginary time oscillation period on energy, by Morais-Smith, Ivlev, and 
Blatter in Ref.EZI where macroscopic quantum tunneling in a dc SQUID has been investigated, by Carriga in Ref.EJ 
concerning the vacuum decai£_in 2 and 3 dimensions in the thin wall limit, by Ferrerat2l for a bubble nucleation in the 
<ft 4 model, and by SkvortsovtLa for the case of the depinning of a massive string from a linear defect. 

It would be interesting to observe a first-orderJ-ransition in an experiment. Measurements of the escape rate of 
Josephson junctions from their zero voltage stateEEl show that the transition from quantum to classical behavior is 
of the second kind. The criterion described above allows us to find metastable systems exhibiting a sharp transition 
without a complicated solution of the equation of motion in the whole temperature range. The problem considered 
in section III C of the depinning of a string with finite mass and friction coefficients is related to the problem of the 



depinning of vortices from columnar defects produced by the irradiation of heavy ions in high-T c superconductors. 
As it has been shown, this problem exhibits a first-order transition. In principle, this sharp behavior can be observed 
in an experiment. However, one should take into account that the point of transition is smeared due to the different 
radii of the defects. Hence, samples with defects of one radius should be taken for an exper,iment. An appropriate 
method to produce defects with identical radii is based on lithographic technique, see Ref. 
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Finally, let us discuss the F — T "phase diagram" of a massive string depinning from a columnar defect. At small 
F, T c increases ~ F and the transition is first-order like. ALJ* 1 < F c , with F c the critical force, the problem exhibits 
a secand-order transition with T c —T ~ y/F c — F, see Ref.Ea. Consequently, there is a "tricritical" point C (see also 
Ref.EJ), where the nature of the transition changes, see Fig. 9. 



T 



classical 




activation 






o 






f \ / quantum 




/ tunneling \ 





F/F c 



1 

Fig. 9 Possible F—T "phase-diagram" of a massive string depinning from a columnar defect. In C the kind of transition changes from 
first- to second-order, i.e., it is a tricritical point. In the vicinity of the point F/F c = 1 the curve is dashed as in the limit F — » F c the 
barrier disappears and the semiclassical approximation is not applicable. 

In the vicinity of the point F/F c = 1 the transition curve is shown to be dashed as in the limit F — > F c the potential 
barrier disappears and the semiclassical approximation is not applicable. 

As it has already been mentioned, the kind of transition is determined by the "competition" of the negative and 

positive eigenvalues of the operator + (1/2) (A — l(2wo)) J , see section [II B . Let us investigate possible spectra 

of this operator for the case of the tunneling of a string. Consider a ID Schrodinger operator h = —ed 2 /dz 2 + Uu n (z) 
of the linearized problem, where the function Uu n (z) is such that the operator h has one negative and one zero 
eigenvalue. Let us try to reconstruct the initial metastable potential V(u). As it has been mentioned in section IIIB, 
the function 'J = duth/dz is a solution of the equation = 0, i.e., uth(z) satisfies the equation 



e ( du 



III 



2\ dz 

where we used the boundary condition u t h, u' th - 
u t h{z) — > z = u7f~(u) we obtain 



-V(u th (z)) = 0, 
if \z\ -» oo and V(0),V"(0) 



(54) 



0. Inverting the equation u 



(55) 



Hence, we have shown that for an arbitrary ID Schrodinger operator of the linearized problem there is only one 
corresponding "original" metastable potential V(u). We have shown therefore that the mapping V(u) — > h is one-to- 
one. Consequently, one can construct a potential V(u) such that the linearized "quantum mechanical" problem has the 
appropriate spectrum. Let the potential V(u) be dependent on some parameter a which might be multicomponcnt, 

in general. Changing this parameter we are modifying the spectrum of the operator [h 1 + (l/2)(/i, — l(2u>o)) J 

and thus we can tune the kind of the transition from quantum to classical behavior. In general the "phase diagram" 
might be more complicated than that in Fig. 9. 

Briefly summarizing, we have studied the behavior of the imaginary time oscillation period of a metastable system 
in the vicinity of the saddle-point solution and have derived a sufficient criterion for a sharp first-order transition 
from quantum to thermal decay of a metastable state. The results have been applied to a comparative study of 
the tunneling of a massive string across a slightly titled periodic potential and the depinning of a string from a 
columnar defect in the presence of arbitrarily strong dissipation. The former problem shows a positive derivative of 
the oscillation period with respect to the amplitude of the motion, in agreement with the second-order transition for 
a purely massive problemu, whereas the latter problem shows a sharp behavior for any mass and friction coefficient. 
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APPENDIX A: CONDITION OF DISAPPEARENCE OF RESONANT TERMS 

Consider the following equation 

ly = hy + s{r,t), (Al) 

where the operators I and h act on time and spatial variables, respectively, s(r, t) is a given vector. The spectrum 
of th e op erator I is negative and the spectrum of h shall be nonnegative except for one eigenvalue. The solution of 



Eq. (Al) takes the form 



-l 



y=[l-h) s(r,i). (A2) 

As [I, h] — 0, the eigenfunctions of the operator (l — h) are products of the eigenfunctions of the operators I and h, 
respectively. The operator (I — h) might have one zero eigenvalue in which case its inverse (I — h) is not properly 



defined. One then can use Eq. (A2) only if the coefficient of the zero eigenvalue eigenfunction of the operator (I — h) 
in the expansion of s(r,i) over the eigenfunctions of the operator (I — h) is equal to zero, i.e., if 

(s,V )=0, (A3) 
where (x, y) is a scalar product and Vq is the eigenfunction of the operator (l — h) corresponding to the zero eigenvalue. 
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